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Abstract 

A variational framework is defined for vertical slice models with three dimensional velocity depending only on 
x and z. The models that result from this framework are Hamiltonian, and have a Kelvin-Noether circulation 
theorem that results in a conserved potential vorticity in the slice geometry. These results are demonstrated 
for the incompressible Euler-Boussinesq equations with a constant temperature gradient in the ^-direction (the 
Eady-Boussinesq model), which is an idealised problem used to study the formation and subsequent evolution of 
weather fronts. We then introduce a new compressible extension of this model. Unlike the incompressible model, 
the compressible model does not produce solutions that are also solutions of the three-dimensional equations, 
but it does reduce to the Eady-Boussinesq model in the low Mach number limit. This means that this new 
model can be used in asymptotic limit error testing for compressible weather models running in a vertical slice 
configuration. 
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1 Introduction 

This paper introduces a variational framework for deriving geophysical fluid dynamics models in a vertical slice 
geometry (i.e. the x-z plane). The work is motivated by the asymptotic limit solutions framework advocated in 
jCu!07j . in which model error in dynamical cores for numerical weather prediction models can be quantified by 
comparing limits of numerical solutions with solutions from semigeostrophic (SG) models. In particular, the SG 
solutions of the Eady frontogenesis problem specified in a vertical slice geometry prove very useful since they can 
be solved in a two-dimensional domain, which means that they can be run quickly on a single workstation. In the 
incompressible hydrostatic and nonhydrostatic cases these solutions are equivalent to exact solutions of the full three 
dimensional equations. As described in |Cul07j . this proves to be a challenging test problem. Using a Lagrangian 
numerical discretisation that utilises the optimal transport formulation, converged numerical integrations of the 
SG model indicate an almost periodic cycle in which fronts form, change shape, and then relax again to a smooth 
solution. However, primitive equation solutions obtained by |GNH92j are rather dissipative due to the need for eddy 
viscosity to stabilise the numerics, and the periodic behaviour is not observed; this leads to a loss of predictability 
after the formation of the front. In |Cul07) . it is suggested that greater predictability in this limit might be possible 
if the numerical solution exhibits energy and potential vorticity conservation over long time periods; it is also 
suggested that a form of Lagrangian averaging may be required to obtain accurate predictions of the subsequent 
front evolution. Since energy conservation can be derived from a variational framework and potential vorticity 
arises from the particle relabelling symmetry, this has motivated us to develop such a framework in the case of 
"slice geometries" in which there are three components of velocity, but they are functions of x and z only. 

Another motivation for our work is that efforts to compare compressible models with the two dimensional SG 
solutions have been thwarted by the fact that it is not possible to construct a compressible vertical slice model 
with solutions that are consistent with the full three dimensional model, with conserved energy and potential 
vorticity. This is because of the nonlinear dependence in the equation of state on the y-dependent component of 
the temperature. Hence, so far asymptotic limit studies of compressible models have only been performed over 
short time intervals corresponding to the initial stages of front formation [Cul08] . In this paper we introduce a new 
compressible slice model that can be used in asymptotic limit studies, since it has a conserved energy and potential 
vorticity. The price to pay is that the solutions are not consistent with the full three dimensional equations. 
However, the model should still be very useful in studying the behaviour of discretisation methods and averaging 
procedures for numerical weather prediction in the presence of fronts. 

Our approach is to derive models in the Euler-Poincare framework [HMR98 . This framework is a way of 
obtaining variational models without resorting to Lagrangian coordinates, by providing formulas that express how 
infinitesimal variations in the Lagrangian flow map correspond to variations in the Eulerian prognostic variables. 
The present paper specialises to the case where all the Eulerian fields are independent of y. This corresponds to 
a subgroup of the group of diffcomorphisms in three dimensions, which can be expressed as a semi-direct product 
of two dimensional diffcomorphisms in the vertical slice and rigid displacements in the y-direction. Having selected 
this group, the Eulcr-Poincarc theory immediately tells us how to perform Hamilton's principle. In this framework, 
the problem of developing slice models reduces to the problem of choosing which Lagrangian to substitute into the 
action. 

The structure of this paper as follows. In Section |2~T1 we identify the slice subgroup, and set up the geometric 
framework. In Section 12.21 we then obtain the general equations of motion corresponding to the Euler-Poincare 
equation with advected density and tracer variables (temperature) . In Section 12.31 we reformulate the equations 
in a more geometric notation, and show that the equations conserve energy in the case of Lagrangians without 
explicit time-dependence; this is shown by recasting the equations in Lie-Poisson form. We also show that the 
equations have a Kelvin-Noether circulation theorem. This circulation theorem differs from the usual circulation 
theorem for baroclinic fluids which have a baroclinic circulation production term on the right-hand side that only 
vanishes if the circulation loop lies on an isentropic surface. In the slice geometry, this baroclinic term can be 
rewritten as the time-derivative of another circulation term, and we obtain conservation of the total circulation on 
arbitrary curves within the slice. This circulation theorem leads to a conserved potential vorticity that turns out 
to correspond to the usual three-dimensional Ertel potential vorticity. We then use this framework to present a 
number of models in the slice geometry. In Section [3] we show how to obtain the Euler-Boussinesq Eady model. We 
present the corresponding Lagrangian-averaged Eady model in Section [4] and introduce our new compressible slice 
model in Section [5j comparing it with the model used in |Cul08 j . Finally we provide a summary and outlook in 
Section [6l The appendices provide proofs and show how this framework relates to known Lie-Poisson formulations 
of supcrfluid models. This relationship is significant since it shows how to build conservative numerical schemes in 
the slice geometry. 
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2 Vertical slice models 
2.1 Definition 

Vertical slice models arise when it is assumed that the forward Lagrangian map takes the form 

4>{X, Y, Z, t) = (x(X, Z, t), y(X, Z, t) + Y, z(X, Z, t)) , (2.1) 
where (X, Y, Z) are Lagrangian labels, (x, y, z) are particle locations and t is time, i.e. 




dY 

Such maps form a subgroup of the diffcomorphismtl]] Diff (ft x R) (where ft E R 2 is the domain in the x-z plane, and 
R represents an infinite line in the y-direction). This subgroup is isomorphic to Diff (Sl)(S)J r (Sl) where (S) denotes the 
scmidirect product, and J~(ft) denotes an appropriate space of smooth functions on f2 that specify the displacement 
of Lagrangian particles in the y-direction at each point in ft. Multiplication in the semidirect product group is 
given by a standard formula, 

(4>u h) ■ (<h, h) = (0i o <t> 2 , 4>i ° h + h)- (2.2) 

The corresponding Lie algebra is isomorphic to X(il)(S)J r (ri) where X(VL) denotes the vector fields on f2, representing 
the two components of the velocity us E X(ft) in the x-z plane, and the smooth function ut £ J~(ft) represents the 
y-component of the velocity. We write elements of X(ri)©J r (0) as (us,ut) where us is the "slice" component in 
the x-z plane, and ut is the "transverse" component in the y direction. In component notation, the Lie bracket for 
the Lie algebra X(ri)(S)J r (f7) of the scmidirect product group Diff (fi)©J r (0) takes the form 

[(us,u T ), {w s ,w t )} = ([us,ws], u s ■ Vw T - w s ■ Vu T ) , (2.3) 

where [us, ws] = us ■ Vws — ws ■ Vw^ is the Lie bracket for the time-dependent vector fields (us, ws) £ 3-(f2), and 
V denotes the gradient in the x-z plane. 

We introduce two types of advected quantities in this framework. 

First, mass is conserved locally, so the mass element Dcfx is advected in three-dimensional space. That is, the 
mass density D(x,y,z,t) satisfies 

(d t + £ {us , Ut) ){D<?x) = (dtD + V • (u s D) + d v (u T D)) <?x = 0, 

with partial time derivative dt — d/dt and partial space derivative d v = d/dy in the y-direction normal to the x-y 
plane. If Ut and D are specified to be y-independent consistently with the slice motion assumption, then the last 
term vanishes and the equation for conservation of mass reduces to advection of an area! density DdS £ A 2 (f2), in 
which D(x, z, t) satisfies the continuity equation, 

dtD + V ■ {u s D) = . (2.4) 

Second, in order to represent potential temperature that has a constant gradient in the y-direction, s = 86/ dy = 
constant, we shall require advected scalars 6(x, y, z, t) that may be decomposed into dynamic and static parts, as 

6(x, y, z, t) = 6 s (x, z, t) + (y - y )s. (2.5) 

Consequently, the three-dimensional scalar tracer equation 

d t s + u s ■ V6> s + u T d y 6 = 

becomes a dynamic equation for 9s(x, z, t) E J~{£i), which satisfies, 

Ms + us- V9 S + u T s = 0, (2.6) 

in which we keep in mind that s is a constant and ut(x, z, t) has been specified to be y- independent. The space of 
advected scalars of this type is isomorphic to J 7 (ft) x R, represented as pairs (9s, s), with infinitesimal Lie algebra 
action 

C(us,u T )(Qs, s) = (u s ■ V6> s + u T s, 0) . 



1 DifEcomorphisms are smooth invertible maps with smooth inverses. 
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2.2 Variational formulation via Hamilton's principle 

in this section we show how to perform variational calculus in the slice geometry. Vector fields of infinitesimal vari- 
ations (ws,wt) in the Lie algebra X(ri)©J r (il) of the scmidirect product group Diff (Sl)(s)J r (Sl) induce infinitesimal 
variations in (us,ut), D, and (6s, s) as follows^ 

5 (u s ,u T ) = (d t ws + [u s ,ws},dtw T + u s ■ Vw T - w s ■ Vu T ) , 

SD = -V ■ (w s D) , (2.7) 
5 (9 S , s) = (-w s ■ V0 S - w T s, 0) . 

For a Lagrangian functional l[(u s ,u T ), (6 s ,s),D) : (X(S)J r (0))(S)(( x M) x A 2 (0)) -» M, we apply Hamilton's 
principle and obtain 

= 5S 

= 5 [ l[(us,u T ),(6s,s),D]dt 
Jo 

si „, A / si p/ „ A / si 



S(us,ut) 

SI \ I SI 

- — , d t w s + (u s ■ V)w s - (w s ■ V)u s ) + { - — , d t w T + u s ■ Vw T - w s ■ Vu T 

OUs I \OUt 

SI \ I SI \ (2 

— , -V • (w s D)) + ( — , -(w s ■ V)(9 S - wts) dt 

-777 V ■ Us ® — - (Vm s ) t - 77— Vu T + £>V— - — Vfl s , w s 

OT OUg \ OM5 J ous Sut oJJ oils 

— V • I u s - — I — — — s, wt ) dt 



f 51 \ 51 \ 



dt Sut 

51 \ I 51 



,Ws) + ( -z , W T 

5us I \ OUT 

where the angle brackets indicate Li inner products with integration over R 2 . The last term makes no contribution 
for velocity variations (ws,wt) that vanish at the endpoints in time. 

Hence, we obtain the Euler-Poincare equations on the slice semidirect product with advected density D and 
scalar 9: 

d 51 „ / 51 \ ,„ xT 51 51 „ nri 51 51 „ 

otdu s V Su s J 5u s 5u T 5D 50 s . . 

d 51 _ / 51 \ 51 [ ' ' 

= r— S . 



dt Sut V fer / 59s 

The system (|2.9p is completed by including the advection equations (|2.4p and (|2.6|) for D and 6*5, respectively. 

2.3 Geometric reformulation and Kelvin-Noether circulation theorem 

Theorem 1 (Energy conservation). If the Lagrangian I has no explicit time-dependence, the energy functional 

h[(m s , tnr), (0s, s),D) = ((m s , m T ), (u s , «r)> - «r), (0s, s), D] . (2.10) 
is conserved for solutions of Equations \2.1$, 12. 6\) and 1(2. 9\) . 



Proof. In Appendix [HI we show that Equations (|2.4[) . (|2.6j) and (|2.9|) arc Hamiltonian, with Hamiltonian given by 
h in Equation (|2.10j) . If I has no explicit time-dependence, then h has no explicit time-dependence and is therefore 
an invariant of the Hamiltonian system. □ 



2 These are standard formulas for defining the variations in Hamilton's principle. Sec HMR98 and Appendix \K\ for details. 
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Theorem 2 (Kelvin- Noether circulation theorem). Equations \2.j$ , i2. 6}) and \2. 9\) imply a conservation law for 
circulation, w, 

in which c(us) is a circuit in the vertical slice moving with velocity us and s = s is a constant parameter 

Proof. The proof of the theorem is facilitated by rewriting the system of equations (|2.4[) . (|2.6j) and (|2.9[) equivalently 
in the following geometric form, 



d \ f 1 61 \ f 1 61 

us 



dt V \D8ut J \D58s 
d 



(2.12) 

. C us J8s + u T s = , 

§- t +£u S ) (DdS) = 0, 

where C us denotes Lie derivative along the vector field us- One may then verify the circulation theorem (12. lip for 
slice models by applying the relation 

4/ v(x,t)-dx=(f ( ^- + C US J (v(x,t) -dx), 

at Jc(u s ) Jc(us)\ at / 

for any vector v(x,t) in the slice. 

Corollary 3. The system of equations &2.12\) implies that the following potential vorticity (PV, denoted as q) is 
conserved along flow lines of the fluid velocity us, 

dtg + us'^Q — for potential vorticity q := — ( curl ( s— - — ) + V#s x V [ — - — (2-13) 

D \ \ D ous J \P> out J J 

Proof. Applying the differential operation d to the first equation in the system (|2.12[) yields 

+ ((curl (~) +i -'V« S Kv(Il)). iiS )=», (2.14) 

where dS is the surface element in the vertical slice, whose normal vector is y. Applying the Lie derivative and 
using the continuity equation for D then yields the local conservation law (|2.13[) . □ 

Upon introducing the new notation, 

1 51 1 51 51 1 51 

Vs: =D5uS' VT:= DS^> *^6D> 1s '-=dW s ' (2 - 15) 



the system (|2.12[) takes a slightly more transparent form 

hr — vt dux — 7s d9s , 

(2.16) 



(d t + Cu s ) ( v s ■ dx) = dn - v T du T - 7s d6 s , 
(d t + C us ) v t = - s 7s , 
(d t + C us )d9 s = - s du T , 
(d t +£ us )(DdS) = 0, 



in which the differential of the third equation has also been taken. Hence, combining the equations in (|2.16|) yields 
the relation 

(d t + C us ) (sv s ■ dx -vt dd s ) - d^ . (2.17) 
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This relation, in turn, yields the Kelvin-Noether circulation theorem as stated above in (|2,11[) . 

^-<f (sv s - v T \78 s ) • dx = <f (d t + £ us )( S v s -dx-v T d9 s ) = <f dn , (2.18) 

* Jc(es) J c(«s) J c(us) 

and potential vorticity conservation as in (|2.13l) . 

dtq + us ■ Vq = for potential vorticity q := — (s curl us + V6>s x Vvt) ■ y ■ (2-19) 

□ 

Remark 4. Note that this circulation theorem is different from the case of general 3D motions, in which the 
circulation is only preserved if the loop integral is restricted to lie on a temperature isosurface. In the special case 
of slice motions, the baroclinic generation term can itself be written as the total derivative of a loop integral. The 
physical interpretation is that q is in fact the usual three-dimensional potential vorticity. Due to the existence of the 
linear y-variation in 9, it is always possible to find an equivalent three-dimensional loop on a temperature isosurface 
that projects onto any given two-dimensional loop in the vertical slice plane. 

Remark 5. In Appendix \B.3\ we will discuss the geometric meaning of the Kelvin-Noether circulation theorem 
\2.11)) and the potential vorticity conservation law V2.13\) from the viewpoint of the Lie-Poisson brackets in the 
Hamiltonian formulation of these equations. 



3 The Euler— Boussinesq Eady model 

3.1 Specialising the Euler Poincare equations to deal with the Eady model 

The Euler-Boussinesq Eady model in a periodic channel of width L and height H , has Lagrangian 

l[u s , u T , D,6,p] = {\u s \ 2 + 4) + Dfu T x + ±.d(z- y \ 9 S + p(l - D) dV, (3.1) 

where g is the acceleration due to gravity, 9q is the reference temperature, / is the Coriolis parameter, and we have 
introduced the Lagrange multiplier p to enforce constant density. We obtain the following variational derivatives of 
this s Lagrangian, 

1 51 1 SI 

vs = t^t — =us, u t = 7^t — =u T + jx, 
L> ous D OUT 

1t = Jd = \ ( |us|2 + U ^ + furx ~ p + t° s { z ~ f ) ' (3 - 2) 

1 61 g ( H\ 51 

Substitution of these variational derivatives into the Euler-Poincare equations in (|2.9p gives 
d t us + u s ■ Vit s + {Vus) T ■ u s + (ut + fx)Vur 

- v (i (ks i 2 + 4) + - J. + i« s (. - f )) - 1 (. - f ) v 6s , m 

d t u T + u s ■ V(it T + fx) = —jj- (^z - y ^ s. 

Upon substituting D = 1, V • us = and combining with equations (|2.4j) and (|2.6p for the advected quantities D 
and 9, the system of equations (|3.3p becomes 

d t u s + u s ■ Vu s - fu T x = -Vp + j-Osz, 

d t u T + u s ■ Vu T + fus ■ x = -j- fz - s > (3.4) 

V-u s = 0, 
d t s + u s ■ V6> s + u T s = 0, 
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where x is the unit normal in the x-direction. 



Remark 6. The system \5.6)) is the standard Euler-Boussinesq Eady slice model. 

3.2 Geometric reformulation and circulation theorem for the Eady model 



Substitution of the variational derivatives in (|3.2[) into the geometric form of the system of Euler-Poincare equations 
in (|2.16l) gives the following equivalent form of this system, 

(d t + Cus) (us ■ dx) = - dp - (u T + fx)du T - 6 S dj s , 
(d t + £u s ) (ut + fx) = -s-ys , , , 

D = l=^V-u s = 0, 

(d t + £ us )dOs = -sdu T . 

Consequently, we recover the Kelvin circulation conservation law (|2.18|) for the Eady model in the form 



d 



f {us - s 1 (u T + fx)V9 s ) -dx= / d(-\us 

Jc(u s ) Jc(us) \ Z 



2 



-P + 7S0S = . (3.6) 



Corollary 7. Equation h3. 6]) and incompressibility imply that potential vorticity (PV, denoted as q) is conserved 
along flow lines of the fluid velocity us in the Eady model, 

dtq + us ■ Vg = for potential vorticity q := (curl us + s^VBs x V(w T + fx)) ■ y . (3.7) 

On denoting u$ = (u,w), Ut = v, this potential vorticity may be written as 

dw du 1 d(v + fx,6') 
q = 1 1 — — -. 

dx dz 9£ Q( x z \ 

dy 



Applying the Legendre transform to the Lagrangian (|3.ip yields the energy 

h[u s , u T , D, 9,p] = jj± (\u s \ 2 + u 2 T ) ~j q D(z- 6 S dV. (3.8) 
Corollary 8. The energy (|3.8j) is conserved for the Eady Boussinesq slice model. 

4 Lagrangian-averaged Boussinesq model 

Numerical forecast models are restricted in grid resolution due to the stringent time requirements of operational 
forecasting, and hence it is necessary to perform some form of averaging on the equations in order to prevent energy 
and enstrophy accumulating at the gridscale, either explicitly by introducing extra terms (i.e. eddy viscosities, 
or Large Eddy Simulation), or implicitly by numerical stabilisation in advection schemes. All of these examples 
amount to some form of Eulerian averaging that leads to dissipation, which is thought to be detrimental to evolution 
of fronts. To avoid this, [Cul07] suggested that some form of Lagrangian averaging may be required, also suggesting 
that it is important for averaging to retain energy and potential vorticity conservation if agreement with the SG 
limiting solution is to be obtained. 

In this section we obtain a Lagrangian averaged Boussinesq model from a variational principle, and so energy 
and potential vorticity conservation will follow immediately. Here, we shall interpret Lagrangian averaging as a 
rcgularisation of the equations that is consistent with the Lagrangian flow map for slice models in Equation (|2.1[) . 
This rcgularisation is obtained by replacing Equation (|3.1[) with 

l[u s , u T , D, 6,p\ = j j (\u s \ 2 + a 2 \Vu s \ 2 + u 2 T + a 2 |Vii T | 2 ) + Dfu T x + j-D (z - j^j 9 S + P (1 - D) dV, (4.1) 
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where a is a regularisation lengthscale. We obtain the following variational derivatives of this Lagrangian, 

1 51 2 2 1 51 o 2 \ 

v s = — - — = (1 - aS/ D )u s , v T = — -z — = (1 - a-V D )u T + fx, 
D dus D out 

31 l /, ,9 9 n g „ ( H\ 

n =SD = 2 C Us| + T ^ + fUTX ~ P+ t 9s V~TJ' (4 - 2) 
1 51 g (H\ SI 



where 



V 2 n = -=-V ■ DV. 



1 



D 

Substitution into the Euler-Poincare equations and applying D = 1 gives 

t 9 
d t us + u s ■ Vu s + Vu s {t s - fu T x = -Vp + 7r#s-2, 

9 t u T + w s • Vu T + fu s -x = (z - — ) s 



V • u s = 0, (4.3) 
d t 6s + us ■ VOs + uts = 0, 

us = (1 - a 2 V 2 )u s , 
ut = (1 — a 2 V 2 )ur- 



This is the Lagrangian averaged Boussinesq Eady slice model. 
Corollary 9. Equations (|4.3[) have conserved energy 



h = J n ^ (M 2 + « 2 |V« S | 2 + 4 + a 2 |Vu T | 2 ) - j-D (z - |) S dV. 
Corollary 10. Equations (|4.3p /lawe Lagrangian potential vorticity conservation 

dtq + us ■ Vg = for potential vorticity q := — (scurlug + W9s x Vut) • y • (4-4) 

5 Sliced Compressible Model (SCM) 

In this section we present a model that is a compressible extension of the Boussinesq Eady model described in the 
previous section. The aim of the model is to provide a framework where nonhydrostatic compressible dynamical 
cores can be benchmarked in a slice geometry. Due to the nonlinear equation of state, it is not possible to write 
down a compressible slice model with solutions that correspond to solutions of the full three dimensional equations, 
and we have to proceed by replacing the full potential temperature 9 in the internal energy by the slice component 
9 S . This approximation would be valid if the potential temperature were slowly varying in the y-direction. However, 
we are more concerned with producing a model that has conserved energy, potential vorticity, supports baroclinic 
instability leading to front formation, so that dynamical cores in this configuration can be compared with the 
corresponding model in the SG limit. 

In the present notation, the Lagrangian for the Sliced Compressible Model (SCM) in Eulerian (x, y, z) coordinates 

is, 

l[u s ,u T ,D,6 s ]= [ ^-(\us\ 2 +u T )+fDuTX + gDz-Dc v 9sUdV, (5.1) 
where LI is the Exner function given by 

, V 

n 



where P q is a reference pressure level and c p and R are gas constants. The equation for an ideal gas becomes 

p IY*- /R = DR9 S U, 
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and differentiating with respect to 9 S and D gives 



<9II d { DRO. 



7-1 



d6 a 86 V Po 

. ,.DRB a ( DR8^~<- 2 
= (7-1)- 



Po V Po 

= 2_l n = Cp 7-n = -^n. 

Similarly wc obtain 

? = — n. 

Note that we use 9 S in both the internal energy term in the Lagrangian, and in the equation of state. This removes 
all y-dependence from the Lagrangian, making a slice model possible. 
We obtain the following variational derivatives of this Lagrangian, 

1 SI 1 SI 

us = — 7 — =us, v T = — 7 — = u T + fx, 
D ous D OUT 

Jp=\ [\ u s? + 4) + f u TX + gz- c p m s , (5.2) 
1 SI _ 1 SI 

where we have used the decomposition (j2.5[) in the last line. 

Substitution of the variational derivatives (|5.2p of the SCM Lagrangian (|5.ip into the Euler-Poincare equations 
in (|2.9[) gives the system 

(d t + C us ) (u s ■ dx) = -c p d s dU + dQ|u s | 2 - gz^j + fu T dx , 

(d t + C us ) ( u t + fx) = sc p U , (5.3) 
(d t + C Us ) s = - su T , 
(dt + Cu s )(DdS)=0. 

Consequently, we recover the expected Kelvin circulation conservation law (|2.18[) for the SCM in the form 

if (u s - s-\u T + fx)VO s ) • dx = . (5.4) 

at Jc(u s ) 

Corollary 11. The system of SCM equations in i5.8\) implies that potential vorticity q is conserved along flow lines 
of the fluid velocity us, 

dtq + us ■ V<7 = with potential vorticity q := — (curl tig + s -1 V0 x V(«t + fx)") ■ y ■ (5-5) 

Corollary 12. These equations are Hamiltonian, with conserved energy 

E = I — (\u s \ 2 + u T ) - gDz + c v DU0' dV. 

Remark 13. The system of SCM equations in \2.11^ may also be written equivalently in standard fluid dynamics 
notation as 

dtus + us ■ Vus — furx = — c p #Vn — gz, 
d t UT + us ■ Vut + fus ■ x = scpU, 
d t 6 s + u s ■ V0 S = — su T , 

(5.6) 

3D 

— + V • (Dug) = • 
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Next we check that the basic state of these equations supports a shear profile (and hence allows barochnic 
instability and frontogenesis). Reverting to more standard notation u s = (u,w), ut = v, the balance equations are 

-fv = -CpO'^, (5.7) 



80 
dy 

dir 



fu = — c p n, (5.8) 



= -c p 6—-g. (5.9) 
oz 

Assuming a x-independent temperature field, then Equation (|5.7[) implies that v = 0. For positive 9, equation (|5.9p 
implies that II will increase with height, and equation (|5.8|l then implies that u decreases with height, leading to a 
shear profile in the basic state. 

We now compare our SCM with the slice compressible model in |Cul08| and identify the differences. On defining 
velocity u = (us, Ut) with us in the vertical slice, and ut transverse to it, the model in |Cul08) in Eulerian (x, y, z) 
coordinates becomes, in the present notation, 

d t u s + u s ■ Vu s - fu T x = - c p 6 VII - gz , 
d t u T + u s ■ Vu T + fu s ■ x = - Cp9U' , 

d t s + u s -V0s = -su T , ( 5 - 10 ) 
dD 

— + \7-(Du s ) = 0. 

Writing the }Cul08| equations in Lie-derivative form yields, cf. equation (|2.17p . 

(d t + C us )(u s ■ dx) = -CpQdU + d^i|w s | 2 - gz^j + fu T dx, 

(d t + C us ){u T + fx) = -Cp9n' , (5.11) 
(d t + C us )9 = -sut, 
(d t + C us )(DdS) = 0. 

These equations differ from the SCM equations in (|2.17l) . by only one term. Namely, the right hand sides of the 
second equation in each set differ, with (— CpII^^) in these equations and (sc p II) in (|2.17|) . It turns out that this 
single difference has important consequences for their respective circulation laws. 

The circulation law for the compressible slice models in |Cul08] is similar to that for the SCM in the previous 
section, but with one important difference. Namely, 

Theorem 14. Circulation for the compressible slice models in |Cul08j is not conserved. Instead, we find 

u 

dt Jc(u s ) x Jc(u s ) 

Proof. The proof uses the first three equations in the system (|5.11|) . The middle two equations yield 

(dt + A^x-s- 1 ^ + fx) do) = i d (s-^n'ofl 2 + 4) + fxdu T . 

Combining this formula with the first equation in the system (|5.1ip then yields the circulation law, (|5.12[) . □ 

Corollary 15. Equation f)5 . 12[) implies that potential vorticity (PV, still denoted as q) is created along flow lines 
of the fluid velocity us, as 

d t q + u s ■ Vg = c p L> _1 Vn x V9 ■ y with PV given by q := D^ 1 (curlw s + s _1 V6» x V(u T + fx)) ■ y . (5.13) 
Proof. Applying Stokes theorem to the circulation equation in (|5 . 1 2[) yields 



u s - s~ 1 (u T + fx)y9) -dx = - (f) c p 9 Vn-dx. (5.12) 



d_ 



(curius + s^ye xy{u T + fx)) -yds = c p ynxye- y ds, (5.14) 

IdS(us) ^ ' JJdS(u s ) 

where y dS is the surface element in the vertical slice, whose normal vector is y. Expanding the time derivative in 
(|5.14[) and applying the Lie derivative relation for D in the last equation of the system (|5.1ip , which is the continuity 
equation for D, then yields the local PV evolution equation in (|5.13[) . □ 
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Remark 16. This is the main difference between the SCM here and in |Cul08| . According to Corollary [TT1 the 
potential vorticity in the SCM is conserved and this conservation is a general property of this class of Euler-Poincare 
equations, as given by Corollary [3j In contrast, according to Corollary 1151 the potential vorticity in the model of 
jCu!08j (when viewed as a slice model) is created whenever the gradients of 8 and II are not aligned. 

6 Summary and Outlook 

In this paper we have shown how to construct variational models for geophysical fluid dynamics problems in a 
vertical slice configuration in which there is motion transverse to the slice, but the velocity field is independent of 
the transverse coordinate. (The vertical slice configuration may be taken as the x-z plane. Then the transverse 
coordinate is y.) Any model developed in this framework has a conserved energy, and corresponding conserved 
potential vorticity. The formulation has a number of interesting geometric features, arising from the semidircct 
product structure of the slice subgroup of the group of three-dimensional diffcomorphisms. Firstly, the formulation 
leads to a Kelvin- Noether circulation theorem in which circulation is preserved on arbitrary loops in the slice, unlike 
the usual circulation theorem in which circulation is only preserved on isentropic surfaces. Secondly, as shown in 
Appendix [B] the equations can always be rewritten in terms of a pair of two dimensional momenta, one comprising 
the x- and z-components of linear momentum, and one formed from the temperature and the y-component of linear 
momentum, plus the density. This formulation involving only two-dimensional momenta and density means that 
potential vorticity conserving numerical schemes for the shallow-water equations can be adapted for vertical slice 
problems. 

This work has led to the development of new model equations: a Lagrangian-averaged form of the Eady model 
of frontogenesis and a new compressible model. We plan to use both of these models to investigate how to improve 
prediction of front evolution, following the agenda set out in [Cul07] . 
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Appendices 

A Euler— Poincare semidirect-product formulation 

The advection equations (|2.4[) - (|2.6p for (9s, s) and D may be rewritten in Lie-derivative notation as 

dt (0s, s) = -C( US!UT ) {0 s ,s) = (-us- V6>s - u T s, 0) , 

d t (D dS) = -C us [D dS) = - d\v(u s D) dS . ( ' ' 

The corresponding infinitesimal variations in (9 s, s) and D, in (|2.7|l induced by the Lie-derivative actions of the Lie 
algebra of vector fields X(ri)©J r (r2) are given by: 

<5 (OS, s) = ~£(wg,w T ) ( d s, s) = (-ws ■ V0 S - w T s, 0) , 
SD dS = -C ws (D dS) = - div(w s D) dS. 1 ' j 

The infinitesimal variations in (us, ut) in (|2.7j) may be expressed in terms of the adjoint action in the Lie algebra 
X(fi)(S)J r (f2) of the semidirect-product group Diff (f2)©J 7 (f2). Namely, 

S (u s , u T ) = (d t w s , d t w T ) - ad (ttSi „ T ) (w s , w T ) ^ 
= (d t w s + [u s , w s ], d t w T + u s ■ Vw T - w s ■ Vw T ) ■ 

For a Lagrangian functional l[(u s , u T ), (6 s ,s),D\ : (X®T(n))®((A°(n) xl)x A 2 (ft)) -> R, one defines Hamilton's 
principle using the L 2 pairing, which is denoted as ( • , • } . Hence, inserting the infinitesimal variational formulas in 
(12.71) for (us,ut), (#s, s) and D yields, in semidirect-product notation, 

= ss 

= 5 f l[(us,u T ),(e s ,s),D}dt 
Jo 

rT 



SI 



d \ 
-, -Q t (ws,w T ) - &d {us , UT) (w s ,w T )) 



o \8(us,u T y dt y 1J " k«sW*""*j/ (A.4) 



51 



'(«; s ,«; r ) (0S, s)\ + ~ div(u>s£ ) )^ dt 



<5(0S,s) 

9 5Z ,„ <5Z SI ,„ , ( SI 

dtS(u s ,u T ) K s ' r; d(u s ,WT) <5(0s,s) V d -° 
5Z 

(tug, W T/ 



S(us,ut) 

In comparison, see equation (|2.8j) for the same Hamilton's principle in vector notation. As before, the last term in 
the previous equation vanishes because (ws,wt) vanishes at the endpoints. The ad* notation in (|A.4[) denotes the 
dual of the ad operation with respect to the L 2 pairing ( • , • ) [HMR98 . Explicitly, the L 2 dual of the ad operation 
is defined by 



(^^st£^y {ws > WT) ) = (w^~y ad{us ' UT)iws,WT) )- 

Likewise, the diamond (o) operation is defined in the present notation by the L 2 pairings, 

\wb) o{0s > s) > {ws ' Wt) ) := (wsb~y - W*> (**•'))' 

( (w <>D, °) ' ■= (w> - C ^ D ) = (w> - div (-^)) • 
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Hence, the last equality of (|A.4[) yields the Euler-Poincare equations on the dual Lie algebra (X(fi)(S)-7 r (^))* with 
the advected areal density D £ A 2 and advected scalars (9 s, s) £ A x 1 in semidirect- product form, as 

d SI .„ 51 51 



ad. 



51 



The system (|A.7[) is completed by including the advection equations (|A.1|) for Z? and (0s, s) 



(A.7) 



B Lie-Poisson Hamiltonian formulation 

B.l Equations on the dual of (X(f)J r (fi))(S)((A (fi) x E) x A 2 (fT)) 
The Legendre transformation to the Hamiltonian is defined by, 

h[(ms, m T ), (9 S , s), D) = ((m s , m T ), (u s , u T )) - l[(u s , u T ), (9 S , s),D]. 
Therefore, we find the variational relations 

51 5h 5h 51 5h 



(ms, m T ) 



, (us, u T ) 



5(u s ,u T )' ' 5(m s ,m T )' 5(6 s ,s) 5(0 s ,s)' 5D 

Consequently, the system (|A.7I) may be written in terms of the Hamiltonian as 



5D 



^(m s ,m T ) 



ad 



*5h/8(ms..m T )( rn Si' m T) 



— (9 S , S - — D, 

5(9 s ,s) 5D 



The advection equations (|A.1[) for (8s, s) and D are then written as 

d 



dt 

d_ 

dt 



(0s, a) 



-C 



8h/ 6(ms ,?tit) 



(9s, s), 



(D, 0) — -Csh/S(m s ,m T )(D, 0) . 



Hence, the entire system (|B.3|) — (|B.4|) may be written in Hamiltonian form as 



d_ 

dt 



(m s ,m T ) 




(0S,s) 




. (AO) . 





adorns, m T ) Do(e s ,a) no(D, 0) 
C a (9 s ,s) 

C a (D, 0) 



5h/5(ms, tut) 
5h/5(9 s ,s) 
5h/5(D, 0) 



(B.l) 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



in which the box □ indicates the appropriate substitutions. The matrix operator in (|B.5j) defines a Lie-Poisson 
bracket dual to the semidirect product action (X(S)J : '(^))(S)((A (r2) x E) x A 2 (fl)) with coordinates (us,ut) € 
X(§)T(n), (9 S , s) G A°(fi) x K and D e A 2 (£l). This identification of the Lie-Poisson bracket with the dual of a Lie 
algebra action guarantees that it satisfies the Jacobi identity. Explicitly, the Lie-Poisson bracket is the following 



{f,h} 



5f/5(m s ,m T ) 
6f/5(9 s ,s) 
Sf /S{D, 0) 



ad* a (m s ,m T ) Do{0 s ,s) □ o (D, 0) 
C a (9 s ,s) 

C a (D, 0) 



5h/5(ms, tut) 
5h/5(9 s ,s) 
5h/5(D, 0) 



(B.6) 



where ( ■ , • ) denotes the L 2 pairing. 

Expanding out the operations in (|B.6j) makes it clear that this Lie-Poisson bracket has the required property of 
being antisymmetric under exchange of / and h. That is, {h, /} = — {/, h}, which is evident upon expanding out 
the operations to express the bracket in (|B.6[) equivalently as 



{/,*} 



(m s ,m T ), 



St 



5h 



5(ms,m T y 5(m s ,m T ) 
5h 



-a 



Sf 



IS, s), £ S f/8{ms,rn T ) 5(Q Si s) ^Sh/5(ms ,m T ) g\ 



(B.7) 



(D, 0), C 



5h 



c- 



Sf 



8f/6(ms,m T )§fT)^ Q) 6h/5(m s ,m T ) Q) 
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Here, C + denotes the L 2 adjoint of the Lie derivative C. In particular, upon denoting 6f/5(ms,rnr) = (u>s,Wt), 
we find the following relations among the operations C + , £ and o, 

/ Sh \ / Sh \ I Sh \ 

((^ s )'4> s ,^)^)/ ■= \w^r C{ws - WT) {dS7S) ) =: \w^j° {9s ' s) ' {ws ' wt) ) ■ (R8) 

Remark 17. If desired, one may now substitute the expressions for Lie derivative (|A.1[) . ad* (|A.5|) and diamond 
(o) (|A.6[) into the L 2 pairings (|B.6[) or (|B.7|) to find the Lie-Poisson bracket {/, h} as an integral over the slice 
domain, O, involving ordinary vector calculus operations. However, the present forms (|B.6[) and (|B.7[) readily reveal 
its semidirect-product nature and suggest further rearrangements, which we pursue next. 



B.2 Equations on the dual of X(S)(A © A 2 © A ) 

To explore the particular case at hand further, one may rewrite the system of equations (|2.4j) . (|2.6j) and (|2.9 
cquivalcntly as, 



d SI 
dt Sus 

d 6i 

dt 5ut 

8 9 





51 

ad - — 

s Sus 



SI 



V9, 



SI 

Sut 



Vii T + DV 



5l_ 
~SD ' 



SI 

Suj 



SI 



50 



s, 



(B.9) 



-C US 9 S - u T s, 



Of 



D = -C US D, 



where C us denotes Lie derivative along the vector field us and we have identified C us and ad* s when acting on the 
1-form density 51 /Sus in the first equation. For more details in this matter, see |HMR98j . 
We define the Legendre transformation to the Hamiltonian in this case by 



h[m Sl m T , 6*5, D; s] = (ms, us) + (m T , u T ) - l[us, u T , ^s, D; s] 



(B.10) 



where the semicolon [ . . . ; s] denotes parametric dependence on the constant s g 
(|B.10j) yields the variational relations 



m s 



51 

5u s 



Sh 



u s 



mT 



51 

Sut 



Sh 



UT 



Sh 



5ms ' Sut ' SrriT 

Consequently, the system (|A.7[) may be written in terms of the Hamiltonian as 



51 

S9s ' 



The Legendre transformation 
(B.ll) 



Sh 
SD 



5J_ 

SD 



dt 
d 



di ms 



= — ad 



Sh 



Sh 



5hf Snis 



Sh 

m s - m T V- + tt-VOs - DW — 

SD 



SmT 58t 



dt 

d_ 
dt 

d_ 

dt 



mT 
Os 
D 



-Sh/5ms m T 



Sh 



-c 



6h/Sm s VS 



Sh 
SmT 



■s , 



s, 



(B.12) 



£-6h,/Sm s D ; 



The corresponding Hamiltonian matrix is 





m s 




ad^ms 


□ mT 


□ o 6 S 


DoD" 




5h/5ms 


d 


mT 




C D m T 





— s 







ShjSmT 


dt 


Os 






s 










5h/5Bs 




D 




_ C n D 













5h/5D 



(B.13) 



in which the box □ indicates the appropriate substitutions. 

After this rearrangement, one recognises (|B.13I) as the Hamiltonian matrix for the Lie-Poisson bracket on the 
dual of the semidirect-product Lie algebra X(S)(A° A 2 © A ) with a symplectic two-cocycle between mj and 9s- 
The Lie bracket for this semidirect-product algebra is 



[(X,f,u,g), (XJ,Q,g)] = ([X,X], X(f)-X(f), X(Q) - X(uj), X(g) - X{g)) , 



(B.14) 
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where, e.g., X(f) = £xf denotes Lie derivative of / by vector field X. The dual coordinates are: ms dual to 
X G X; niT to / G A°; 9 s to ui G A 2 ; and D to g G A°. The spaces in which the coordinates themselves are defined 
are (ms, mr, &s, D) G (A 1 ® A 2 , A 2 , A , A 2 ) and s G M is a parameter. The second part of the bracket (|B.13[) is the 
standard two-cocycle (symplectic form) on A°ffiA 2 arising from the natural projection X(S)(A ©A 2 ©A°) — >• A°©A 2 . 

Remark 18. The Hamiltonian matrix with the two-cocycle in (|B.13[) has been seen before. Namely, it is the same 
as that for 4 He superfluids [DV80, HK82] in the spatially two-dimensional case. For 4 He superfluids, the function 
9s here plays the role of the phase of the Bose-condensate wave function, whose gradient V#s is the superfiuid 
velocity. The other variables ms, mr and D correspond respectively, to total momentum density, mass density and 
entropy density of the superfiuid. 



B.3 Equations on the dual of Xi©(X 2 © A ) 

[HK82] showed that the two-cycle in (|B~13l) may be removed by transforming to new variables 

(ms,mT,9s,D) ^ (ms 7 mR,D) where m R := (s)~ 1 ititV9s ■ 



(B.15) 



The quantity mji is the momentum map for right action of the diffcomorphisms on the buoyancy 9s in two spatial 
dimensions, see |HM04j for more details. The resulting Lic-Poisson bracket has the standard form dual to the Lie 
algebra X\(§)(X 2 © A ), whose Lie bracket is 



[(X 1 ,X 2 J), (X l7 X 2 ,f)] = 

([Xi.Xr], [X 2 ,X 2 ] + [X U X 2 ] - [X U X 2 ], X^f) 



(B.16) 



Dual coordinates in this case are: ms dual to X\ G 3£i; tor to X 2 G 3L 2 ; and D to / G A . 

Transformation of the Hamiltonian matrix (|B.13[) into these variables yields the following Lie-Poisson Hamilto- 
nian system 



d_ 

di 



m s 




m R 




D 





adpTOs ad^mjf Do D 
adpTOii ad^mij 
C n D 





Sh/5ms — 


u s 




Sh/SniR = 






Sh/SD = 


P 



(B.17) 



This system produces a system of equations for relative momentum (mg — m^), momentum map — (s) 1 mT^9s 
and mass density D, given by 



d t (m s - m R ) = - ad* s (ms ~ m R ) -poD , 
d t m R = -ad\ Us+UR) m R , 
d t D = -C us D. 



(B.18) 



Upon evaluating p o D = DVp, the first of these equations explains the geometric origin of the Kelvin-Noether 
circulation theorem (|2.11| that was found by direct manipulation in Section 12.31 Together, the three equations in 
(|B.18j) show that the slice dynamics may be expressed in terms of (ms, m R , D) as a Lie-Poisson Hamiltonian system 
on the semidirect product 

DifFi(fi)©(Diff 2 (fi) xA 2 (!!)), 
in the slice domain Q. When D = 1 is imposed, we have V • us = and this simplifies to 



SDiffi(n)©Diff 2 (n) 



